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Abstract: We provide new expressions for the scattering amplitudes in the soliton- 
antisoliton sector of the elliptic sine-Gordon model in terms of cosets of the afHne Weyl 
group corresponding to infinite products of q-deformed gamma functions. When relaxing 
the usual restriction on the coupling constants, the model contains additional bound states 
which admit an interpretation as breathers. These breather bound states are unavoidably 
accompanied by Tachyons. We compute the complete S-matrix describing the scattering 
of the breathers amonst themselves and with the soliton-antisoliton sector. We carry out 
various reductions of the model, one of them leading to a new type of theory, namely an 
elliptic version of the minimal D^Xj-affine Toda field theory. 



1. Introduction 

By investigating a Z4-symmetry of the particle wavefunctions for a soliton and an antisoli- 
ton, the elliptic sine-Gordon model was introduced originally by A.B. Zamolodchikov more 
than twenty years ago |IJ. The S-matrix was found to correspond to the transfer matrix 
of Baxter's eight-vertex model [§, ||, |J]. The model has two free parameters v and £ which 
were mutually restricted in Q in order to avoid the presence of Tachyons. Thus, if it was 
not for that restriction, the model could be viewed as a generalization of the sine-Gordon 
model. Unfortunately, when constraining mutually the parameters, it corresponds in the 
trigonometric limit only to the soliton-antisoliton sector with no bound states. This means 
that the entire breather sector is absent. 

The purpose of this paper is to investigate whether a meaningful breather sector of the 
model can be constructed. When relaxing the constraint on the parameters several poles 
in the soliton-antisoliton S-matrix amplitudes will move into the physical sheet. We will 
demonstrate below that some of them lie on the imaginary axis and therefore, with the 
appropriate sign of the residue, admit the usual interpretation as bound states correspond- 
ing to breathers. In addition, there are redundant poles which are of a tachyonic nature 
as they are in the physical sheet beyond the imaginary axis. In the context of integrable 
models solutions for scattering amplitudes with such properties have been discarded right 
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away up to now. Nonetheless, Tachyons have emerged initially undesired in other areas 
and subsequently have been turned into virtues. For instance, recently there has been 
great activity in the context of string theory || ||, [/], ||, ||, [[(], |TT], [l^, [U], [l5|, |0| , where 
Tachyon condensates have turned out to be important. They also appear to be very useful 
in the context of cosmological considerations [17, 18, 19, 20, 21, 22, 23, 24, 25]. Here we 
want to adopt the lesson from string theory and cosmology and allow them to be present. 
The gain of this attitude is that we have additional poles at our disposal which admit the 
usual interpretation as bound states, which can be associated to breathers. 

Our manuscript is organized as follows: In section 2 we assemble and derive various 
properties of q-deformed gamma functions and relate them to Jacobian elliptic functions. In 
section 3 we demonstrate how the affine Weyl group can be utilized to solve the central func- 
tional equations leading directly to an infinite product solution for the soliton-antisoliton 
backward scattering amplitude in terms of q-deformed gamma functions. In section 4 we 
discuss the soliton-antisoliton sector. Using a slightly modified procedure proposed origi- 
nally by Karowski and Thun [26, 27], the soliton-breather amplitudes are constructed in 
section 5 and the breather sector is discussed in section 6. By specifying the parameters 
to certain values we reduce the model in section 7 to several models. Our conclusions are 
stated in section 8. 



2. q-deformed gamma functions and Jacobian elliptic functions 

Many of the well studied integrable quantum field theories which allow for backscattering, 
such as the sine-Gordon model, are known to possess a generic form for their scattering 
amplitudes which consists of infinite products of Euler's gamma functions 

OO P r(Q I T i \ P P 

fc=0i=l v y ' i=l i=l 

The rapidity dependent prefactor f{9) is usually some finite product of ratios of trigono- 
metric functions and the values of Xi,yi,a,p are specific to each individual model. The 
constraint on the sums of the Xi , yi is a necessary condition for the convergence of the 
infinite product (a proof of this can be found e.g. in |28|| ). Below we will describe how the 



structure (2A) can be generalized very naturally by replacing in a controlled way the usual 
gamma functions by their q-deformed counterparts and the trigonometric functions in the 
prefactor by their elliptic version. 

q-deformed quantities have turned out to be very useful objects as they allow for 
instance to carry out elegantly (semi)-classical limits when the deformation parameter 
is associated to Planck's constant. In the elliptic sine-Gordon model we have the two 
free parameters v £ [0, oo) and I G [0, 1] at our disposal. The former is the analogue 
of the coupling constant in the sine-Gordon model 1 and the latter is the modulus of the 
Jacobian elliptic functions. It is the 0. which is associated most naturally to a deformation. 



1 This motivates to take v to be positive, as negative values of it correspond to a regime in which according 
to the arguments of Coleman |E9| the ground state of the theory is not bounded from below. 
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Accordingly, we define a deformation parameter q and its Jacobian imaginary transformed 
version, i.e. r — * — 1/r, as 

q = exp(z7rr), g = exp(— iir/r), T = iK\^pjKt . (2.2) 

We introduced here the quarter periods Kg of the Jacobian elliptic functions depending on 
the parameter t £ [0, 1], defined in the usual way through the complete elliptic integrals 
Kg = Jo 2 (1 — ^sin 2 6)~ l l 2 d6 . Recalling the well-known properties 

lim Kp = lim K\_p = ir/2 and lim K.\—i = lim Kp — ► oo , (2.3) 

e-*o t-*x i-+o t->i 



the definitions fl2.2| ) obviously mean that the trigonometric limits correspond to the "clas- 
sical" limit in the variables q, q as 

lim = lim and lim = lim . (2-4) 

£-►0 q-fl t-*\ g-»l 

It will turn out below that quantities in q will be most relevant for our purposes and 
therefore we state several identities directly in q, rather than q, even when they hold for 
generic values. The most basic q-deformed objects one defines are q-deformed integers 
(numbers), for which we take the convention 

(2.5) 



q- 1 



They have the obvious properties 



limlnU = n, (2.6) 

lim ^ + — I 1 for m, m' 7^ ^. 

l— >o [n' + m'r]g 1 n/n' for m = m' = 

With the motivation in mind mentioned at the beginning of this section we define a q- 
deformed version of Euler's gamma function 



-n [x+n] . [n] , • 



(2.8) 



The crucial property of the function Tg, which coins also its name, is 

UmT„(x + 1) = MmTs(x + 1) = TT — ^— (^-^] = T(x + 1) . (2.9) 

n=l 7 

We report now various properties of this function which will be useful below. We can relate 
deformations in q and q through 

g -(*+T/2-i/2)2 T q (y)T q (l-y) _ T q{ -y/ T )T q (l + y/T ) 
$(»+T/2-l/2) a r$(x)r $ (l-s) T q (-x/T)T q (l + x/r) ■ { - > 
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Frequently we have to shift the argument by integer values 



r^x + i) = ^- i [x]^r 



(2.11) 



Relation ( 2.11| ) can be obtained directly from fl2.8| ). As a consequence of this we also have 



m— 1 



T q (x + m) = T q (x) 11 q x+l - l [x + I] 
1=0 

m—1 

T q {x) = T q (x-m) Yl f 



zx-l-2 



X-l- l]g 



m G Z 



m G 



1=0 



Whereas (2.11)-( |2.12[ ) hold for generic q, the following identities are only valid for q 



r$(l/2-T/2)r$(l/2 + r/2 
T q (x + 2t 



n 



T q (y + 2r 
T q {Xi)T q {Xi ± t/2 



r q (xi ± t/2 



lim I 



^ 1/4 r^(i/2) 

Tg(x) 

rr Fg 2 (xj) 



(2.12) 
(2.13) 

i 

(2.14) 
(2.15) 



lim 



1 



7T 



for 



i=l i=l 

P V 

i=l i=l 

x / (2.18) 



The constraint on the sums in ( |2.16| ) and ( p. 17 ) was already encountered in fl2.1| ). Most 
of these properties can be checked directly by means of the defining relation ( |2.8[ ). We 
comment on some of the derivations below. 

The singularity structure will be important for the physical applications. It follows 
from (|2.8j ) that the T^-function has no zeros, but poles 



lim 

9 F>=« 



r$(0 + l) ->oo 



for m G Z, n £ N . 



(2.19) 



Furthermore, we will employ the Jacobian elliptic functions, to generalize the prefactor f(9) 
in ( |2.lD , for which we use the common notation pq(z) with p,q G {s,c,d,n} (see e.g. [Q 
for standard properties). We derive important relations between the q-deformed gamma 
functions and the Jacobian elliptic sn-function 



sn (x) 



i ^(jgj =f S)r*(i 



x _i_ r N 
2K* 31 2/ 



r <?(2K (! ) r <?( 1 2^) 

1 i ix \ "P / 1 
) L 9V2 



21^ _ 



r,(i 



2-RTi 



;)r 9 ( 



(2.20) 
(2.21) 



IK, 



These relations can be used to obtain some of the above mentioned expressions. For 
instance, recalling that sn(i'Q) = 1, we obtain ( |2.14| ). With (2.8) we recover from this the 
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well known identity sn(iKi^/2) = i/i 1 ^ . The trigonometric limits 
lim sn(x) 

lim sn(x) 
l->\ 

can be read off directly recalling ( |2.3[ ), ( |2.9[ ) and presuming that ( |2.18| ) holds. We recall 
the zeros and poles of the Jacobian elliptic sn(#)-function, which in our conventions are 
located at 

zeros: 0*£ o = 2lK e + i2mK 1 ^ i l,m€Z (2.24) 

poles: 6 l ™ p = 2lK e + i(2m + l)K x _ t l,meZ. (2.25) 

We have now assembled the main properties of the q-deformed functions which we shall 
use below. 



lim sn(x) 
<2->l 



lim sn(x) 



7T 



sin(x) 



r(2)r(i - 2" 
i r(i + f)r(i-f) 
i r(i-f)r(f) 



tanh(x). 



(2.22) 
(2.23) 



3. AfRne Weyl group and the unitarity /crossing relations 

The functional relations of crossing, unitarity and bootstrap for the scattering amplitudes 
are usually solved by means of Fourier transformations, thus leading in general directly to 
integral representations for the desired quantities. When solving these integrals one ends 
up with infinite products over gamma functions of the type ( f2.1| ) for scattering amplitudes 
in non-diagonal theories or trigonometric functions when backscattering is absent. For 
the elliptic sine-Gordon model so far only the analogue of the integral representation was 
presented [l[ in form of a discretized version. Instead of solving the discrete integrals in this 
function we provide here a systematic procedure which leads directly to product solutions 
by utilizing the affine Weyl group W(g). 

Let us first assemble the necessary mathematical tools and jargon. In general an affine 



Weyl reflection related to a simple root ai, of a Lie algebra g [31| may be realized by the 
map 

(?i,n (x) = x - (x ■ ai + na-i , (3.1) 

where hi denotes a coroot and n an arbitrary integer. In particular for n = one recovers 
the ordinary Weyl group W(g) and for x = the translation group on the coroot lattice. 
Hence the affine Weyl group may be thought of as a direct product 

W(g) = W(g)®T , (3.2) 
with T denoting translations on the coroot lattice. For the case g = A\ (W (g) ~ the 



infinite dihedral group) one has only one simple root and (3J.) becomes 

a n (x) = na — x . (3-3) 
This group may be generated by two generators 

o\ (x) = a — x and o"o = — x , (3-4) 
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having obviously the property a 2 = <7q = 1. Defining then the transformation 

a := <7i<7o (3-5) 

one has 

a n {x) = x + na . (3.6) 

We denote by Nq,N\ the number of times the generators ctq,cfi occur in an arbitrary 
element of W (A\). Then two types of subgroups {a 2n , N Q ] , {a 2n , Ni] C W (A\) are 
constituted by the elements with Nq and JVi even, respectively. The right cosets of 
{a 2n , Ni} may be divided into those with Ni even {a 2n I, N\] , {a 2n ao, Ni} and Ni odd 
{<7 2n+1 J, A^i} , {<7 2n <7i, iVi}. Similarly one may divide the right cosets of {a 2n ,No} into 
those with N even {a 2n I,N } , {cr 2n er i, iVo} and N odd {a 2n+1 I,N } , {a 2n+1 a!, N } . 

Assuming now that W (A\ ) acts in the complex rapidity plane one may specify a and 
define the "unitarity" and "crossing" transformation on an arbitrary function / (9) by 

*of (0) = f (-0) , o- x f (9) = f (in - 6) . (3.7) 

Then one obtains 

a n f (9) = f(9 + niir) . (3.8) 

We have now provided all the tools to solve the key equations in this context. In Q 
two functional relations for the soliton-antisoliton transmission amplitude (equations (2.10) 
and (3.8) therein), which we denote by c(9), were derived from crossing, unitarity and the 
Yang-Baxter equations. Once this amplitude is known, some simple relations provided in 
Q suffice to construct the remaining ones in the soliton-antisoliton sector. The equations 
to be solved are 

c(9) = c(m - 9) (3.9) 

<0)c(-9) = 2, TT^4rr/ \ ■ ( 3 - 10 ) 
sn z yK/v) — sn z (it) jv) 

We make now an ansatz by taking the ratio of right cosets in which Nq is even and odd 
r(e] _, {^ n I,N }{a 2 -a 1 ,N } " a 2 * p (9) a^p (9) 

C W K { a 2n+l I}No y{ a 2n+l aiiNo } ^ ^ 2fc+1 p (9) <T^G X p (0) ' K ' 



oo 



n p[9 + 2mk\p[-9 + 2iti(k + 1/2)] 
K l} 1 p[e + 2TTi(k + l/2)]p[-9 + 2iii(k + l)\ ' { ' 

Here k is a constant and p (9) an arbitrary function which remains to be fixed. One observes 
that the crossing relation fl3.9| ) is solved by construction, whereas ( 3.10| ) requires that 



hc p(9 + 2m)p(2m 



rg[-j]r|[l + jjTgg - - 9 + ^[-9 - |]Tg[l + 9 + |] 

r i[-|] r l[i + i W - 5]rg[i - § + §]r$[-0 - §]i>[i + + §] 



We replaced here in ( |3.10| ) the sn- by I^-functions using some standard identities for 
Jacobian elliptic functions together with ( p.20| ) and abbreviated for compactness 

A = -n/K e v and 9 = i9/2K e u. (3.13) 
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The problem of solving (|3.9| ) and ( 3,1C| ) has now been reduced to the much simpler task of 
fixing p (9) and k. We find 

Eghj]ij[i+j] p m (0) 

fc = t t - and p.- AO) = j-r , (3.14) 

where p n> i(9) and Pdj(d) could be any of the functions 

p n>1 (0 + 2™) = r 4 [e - |]r 4 [i - + 1] ^(0 + 2™) = r$[0 - |]r^[i - + §], 
Pni2 (0 + 2™) = r# - |]r 4 [i + + 1] Pd>2 (0 + 2™) = r$[$ - |]r^[i + £ + §], 
Pn , 3 (0 + 2™) = r 9 [-e - |]r*[i + 3 + |] ^ + 2™) = ^[-0 - §]r^[i + + §], 

Pni4 (0 + 2tu) = T 4 [-8 - - + |] ^ 4 (0 + 2tu) = r 4 [-e - Z]T 9 [1 - + |]. 



Hence, we notice that the solution of (^^) and ( |3,10| ) is by no means unique and there exist 
additional ones to the one presented in jlj] in form of a discrete integral representation. The 
function in the numerator p ni (0) may be selected by the requirement that we would like 
to obtain the corresponding quantities of the sine-Gordon model in the trigonometric limit 
and Pd,j(6) by discarding solutions which yield undesired poles. We end up with the choice 
Pa 4(0) = PnA{®)/Pdi{®)i which after substitution into ( |3.12 ) yields. 



(e) = Eg^hjM-f] yr r 4 [e - kX\r^[i + e-(k- i)A] 



r 9 [i + $]t 9 [-§] j£ r 4 [-e - k\]r 9 [i -e-(k- i)x] 
rrf-e - ( k - |M]£g[i -e-(k- |)A]r g [i + 9-kx+ T 1 ] 

T . [e — (k+ i)A]r 5 [l + - (fc - |)A]r 5 [l -9 - (k - 1)A + 5] 

r# - k\- r\r 4 [-e — (k— i)A - §]r^[i - - (& - |)a + §] 



(3.15) 



^[-0 - (k - i)A - |]r*[0 - (k - |)a - |]r 9 -[i + - (*, - |)a + 1] 

We want to conclude this section with a general remark on the method provided to solve 



the functional relations ( |3.9[) and ( 3.1C ). Of course we could have started right away with 
the ansatz ( 3.12 ) instead of introducing the affine Weyl group in the first place. However, 
this formulation automatically supplies a certain systematic. To illustrate this further we 
provide another example of some important functional relations occurring in the context 



of integrable models, that is Watson's equations, see e.g. |32|| , for the minimal form factors 

F min (9 + m) = F min (ivr - 9) and F min (9) = F min (-9)S(9) . (3.16) 

Making here an ansatz by taking the ratio of right cosets in which JVi is even and odd 

_ {a^I, Nl } KVVj _ fr a^p(9)a^a 0P (9) 
wW ~ K {a^I,N l }{a^a 1 ,N 1 }~ K j^a^- 1 p(0)a^a 1 p(dy { ) 

we see that the first equation in (|3.16[) is solved by construction whereas the second requires 
that S(9) = p (9) o\~p (9) /<JoP (9) ap (9). Having a concrete expression for S(9) one can now 
easily determine p and p. For more complicated functional relations one may use groups 
of higher rank. 
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4. Soliton-antisoliton sector 

Let us now discuss in more detail the elliptic sine-Gordon model. Scattering amplitudes 
are obtainable in general from the computation of matrix elements, but it is also well 
established that in an integrable (=factorizable) theory they may be derived equivalently 
by analyzing the Zamolodchikov algebra. Its associativity corresponds to the Yang-Baxter 
equations. Its generators are thought of as particle creation operators and therefore internal 
symmetries of the model are respected by this algebra. Considering a theory with two 
particles which are conjugate to each other, say a soliton Z and an antisoliton Z, one 
may demand a ^-symmetry, that is one requires invariance under Z — > exp(z7r/2)Z, Z — > 
exp(— m/2)Z. The most general version of the Zamolodchikov algebra respecting this 
symmetry then reads [j]J 

z{e l )z(9 2 ) = a{e l2 )z(e 2 )z(9 l ) + d{9 12 )z(e 2 )z(e 1 ) , (41) 
z(0i)z(0 2 ) = b(e u )z(e 2 )z(ei) + c{e l2 )z{e 2 )z{e{) , (4.2) 

with rapidity difference 6\ 2 = 6\ — 2 . The charge conjugated relations also hold, that is 
Z «-> Z. In comparison with the more extensively studied sine-Gordon model the difference 



is the occurrence of the amplitude d in (4.1), i.e. the possibility that two solitons change 



into two antisolitons and vice versa. Invoking the Yang-Baxter equations, crossing and 
unitarity one finds the following solutions for the amplitudes 

,„ = m f\ i>[£ - (k + i)A]i> [i + e - k x]T, P [-e-i^\\T f [i - e - a^x\ 

d(0) = -Vlsn{i6/v)sn(Tr/v)a(e) = d(0)a(6), (4.6) 



$(0) = 91 \ 2i qy 21 5 [ 2 \ 2J gL — • (4-7) 



Here we use as a common factor a(9) rather than c(0) for which we have explained in the 
previous section how to solve the key functional equations ( |3.9|) and ( |3.10|) . We also used 
relation ( 2.11[ ) to simplify ( |3.15 ) and refer the reader to [|| for the details on how to relate 



the different amplitudes to each other. 

Our solutions (4.3)-(fO|) differ in form from the one presented in M, where a discretized 



version of an integral representation was provided. Instead we derived here directly an 
infinite product representation in the previous section, which is a natural generalization of 
a very common version, of the form (|2.l| ), used in the context of the sine-Gordon model. 
One of the advantages of our formulation is that it exhibits very explicitly the singularity 
structure. Furthermore, it is easier to handle under shifts of the argument as in integral 
representations such shifts will often be prohibited by convergence requirements. 
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4.1 Singularity structure 

In order to extract the singularity structure for the amplitudes ( |4.3| )-( |4~6| ) we recall the rela- 
tions fl2.19|) , ( |2.24[ ) and ( |2.25| ), which suffice to read off the poles and zeros of a(9),b(9), c{6) 
and d(9). We depict them most conveniently in a figure. 
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Figure 1: Poles of the soliton-soliton amplitude a(9) = b(iir — 9) for the concrete values v = 0.25, 
I — 0.344, that is vK\-i ~ 0.504 and vKi ~ 0.435. The arrows indicate in which directions the 
"strings" of poles are extended for growing integer values n or m. 

Explicitly, we just report here the poles which have potentially a chance to be situated 
inside the physical sheet, i.e. < Im6* < n. We take !,raeZ,n6fJ and associate always 
two sets of poles 6™ p and #™ m p to a(9), 9%™ and 9%™ to 6(0) etc. 

0™ = 2mvK x _ t + i2nuK e , 9 n ™ p = (2m + l)vK X -i + i(vr - 2nuK e ), 

9 l ^ p = 2mvK 1 ^ + i{K-2lvK l ), 9 l ^ p = (2m + l)vK x ^ + i2luK t , 
9 l ™ p = 2mvK x . t + i2luK e , 9 l ™ p - ' ' ' ~ ' ' 1 



2mvK x _ t + i2luK e , 9 l ™ p = 2muK 1 _ e + i(ir - 2luK t 

,™ p = (2m + V)vK x _ t + i2luK £ , 9 l ™ p = (2m + l)uK x _ e + i(n - 2nrh, 



One readily sees from ( ^4.8| ) that if one restricts the parameter v > tt/2K£ all poles move out 
of the physical sheet into the non-physical one, where they can be interpreted in principle 
as resonances, i.e. unstable particles. This was already stated in [Q| where the choice 
v > ir/2Ki was made in order to avoid the occurrence of tachyonic states. In fact, even 
in that regime there are Tachyons present, since the poles in the non-physical sheet with 
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negative real part can not be explained on the basis of the Breit-Wigner formula. Their 
occurence can be avoided by an additional breaking of parity (see discussion in ]33| ) . The 
restriction on the parameters makes the model somewhat unattractive as this limitation 
eliminates the analogue of the entire breather sector which is present in the sine-Gordon 
model, such also that in the trigonometric limit one only obtains the soliton-antisoliton 
sector of that model, instead of a theory with a richer particle content. For this reason, 
the arguments outlined in the introduction and the fact that the constraint does not yield 
any Tachyon free theory anyhow, we relax here the restriction on v. The the poles 

e h, P = e c2,P for < n < n ^ax = [ir/2vK e ],n e N (4.9) 

are located on the imaginary axis in the physical sheet and are therefore candidates for the 
analogue of the n^-breather bound states in the sine-Gordon model. We indicate here the 
integer part of x by [x\. In other words, there are at most n max — 1 breathers for fixed v 
and 1. The price one pays for the occurrence of these new particles in the Z4-model is that 
one unavoidably also introduces additional Tachyons into the model as the poles always 
emerge in "strings". It remains to be established whether the poles ( fl.9| ) may really be 
associated to a breather type behaviour. 



5. Soliton-breather amplitudes 



Once the solitonic sector of a theory is constructed there is a well defined bootstrap proce- 
dure proposed by Karowski and Thun 



27 1, which allows to complete the theory and to 
construct also the breather sector. There are numerous solutions known for a non-trivial 
soliton sector, since for instance for all affine Toda field theories with purely imaginary 
coupling constant the Yang-Baxter equations can be solved by representations of the cor- 

Nonetheless, the completion of the models has 



responding quantum group |34|, £5 



not been carried out for very many cases [37, [3^, 39]. In reverse, one should note that 
there exist many scalar theories which do not possess a known solitonic counterpart. One 
reason for that is, that the procedure [26, ^7]] is not reversible and one can in general not 
construct the solitonic sector from the breather sector alone. Drawing a loose analogy to 
group theory one may think of the breather sector as a subgroup, which of course does 
not contain the information of a larger group in which it might be embedded. It is very 
desirable to complete the picture, tie up the loose ends and construct the respective missing 
sectors. 

To start the construction, one first of all has to find creation operators for the parti- 
cles corresponding to the bound state poles in the soliton sector. We presume here that 



the particles related to the poles (4.9) are breathers and borrow some intuition from the 
classical theory, exploiting the fact that a breather is an oscillatory object made out of a 
superposition of a soliton and an antisoliton. For the sine-Gordon model this prescription 



was used in [27] to define the n th -breather particle creation operator. Even though we 
do not have a classical counterpart for the elliptic sine-Gordon model, we follow the same 
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approach here and define the auxiliary state 

z n (9 1 ,e 2 ) -.= i= [z{e 1 )z{e 2 ) + (-i) n z(e l )z(e 2 )] . (5.1) 

Obviously, this state has properties of the classical sine-Gordon breather being chargeless 



and having parity (—1)™. Choosing thereafter the rapidities such that the state (|5.1|) is 
on-shell, we can speak of a breather bound state 



lim Z n {9 1: 9 2 ) = lim Z n {9 l ,9 2 ) = Z n {9) . (5.2) 



Here is the fusing angle related to the poles in the soliton-antisoliton scattering ampli- 
tudes. To simplify notation we abbreviate in what follows a^- = a(9{j), bij = b(9ij), etc. 



With the help of the braiding relations (4J) and ( [4.2[) we compute the scattering amplitude 



between a soliton and the auxiliary breather states (|5.l| ) 

Z n {9 l ,9 2 )Z(9 3 ) = i= [a l3 b 23 Z(9 3 )Z{9 l )Z{9 2 ) + (-l) n b 13 a 23 Z(9 3 )Z(9 1 )Z(9 2 ) 

+c l3 c 23 Z(9 3 )Z{9 l )Z{9 2 ) + (-l) n d l3 d 23 Z{9 3 )Z(9 l )Z(9 2 ) 
+b l3 c 23 Z{9 3 )Z{9 l )Z{9 2 ) + (-l) n c l3 a 23 Z(9 3 )Z(9 l )Z(9 2 ) (5.3) 
+ d 13 b 23 Z(9 3 )Z(9 1 )Z(9 2 ) + (-l) n a 13 d 23 Z(9 3 )Z(9 1 )Z{9 2 )] . 

Going now on-shell in the sense as defined in ( |5.2|) , splitting the fusing angle into 9\ = 
9q— 9 2 i/2, 9 2 = 6*o+ 9 2 i/2 and denoting 9q 3 = the following identities can be shown to 
hold 

h 3 c 23 + (-l) n c 13 a 23 = for 9 n 2X = i(n - 2nvK e ) (5.4) 

di 3 b 23 + (-l) n ai 3 d 23 = for 9 n 2X = i(vr - 2rwK t ) (5.5) 

h 3 a 23 + (-l) n c 13 c 2 3 = a 13 b 23 + {-l) n d 13 d 23 for 9 n 2X = i(n - 2nvK t ) . (5.6) 

Therefore the braiding relation ( |5,3| ) reduces to a diagonal scattering process, i.e. there is 
no backscattering amplitude. We obtain after some algebra 

Z n (9 )Z(9 3 ) = S bnS (9)Z(9 3 )Z n (9 ) , (5.7) 

with 

S bnS (B) = a 13 b 23 + (-l) n d 13 d 23 (5.8) 

sn(i9/i> — tt/2v + nKA r/) 9 ir 9 ( i9 ir _ \ . _ 

[fen 2 -sn 2 — + — + nKi - l]ai 3 a 23 (5.9) 



sn(i9/v + 7r/2^ + nif|) f \ z/ 2f 
where 

r> [i + 9 + 1 - §]r> [-0 - f - §]r> + | + 1 ]r> [e + 



ai 3 a 23 = $i 3 ^23 

A ;„ \ i r £) i « A ;„ \ i \o i n i A r„\i r 2 n A 



i> [i - + 1 - f ]r> §]i> [0 + f + f ]r> [-0 + £ - §] 



X 



[g _ % + a _ fc^-fl + I - I - fc%[0 + f + f - fc%[-0 - | - | - fc% 

/, :0 M-i + r + 1 - 1 - fc%[-0 + f + 1 - fc%[§ - f - \ - k\] e 



n 



x n f\ [9-^ 2 +\-k\ + l\l 2 [-9 + ^-\-k\-l]l 2 
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Similarly, we compute the braiding of Z n (9\,92)Z{9^) which leads us to 

S bnS {0) = S bnS (0) . (5.11) 



Denoting the anti-particle always by an overbar, we deduce from ( |5.11 ) that the breathers 
are self-conjugate due to the general relation S a b = S^, that is b n = b n . 

The matrix Sb nS (0) contains various types of poles, i) simple and double poles inside 
the physical sheet beyond the imaginary axis, which are redundant from our point of 
view as they are of a tachyonic nature, ii) double poles on the imaginary axis which can 
be explained by the usual "box diagrams" corresponding to the Coleman-Thun mechanism 
fiofl , iii) simple poles in the non-physical sheet which can be interpreted as unstable particles 
(see e.g. |33| and references therein) and iv) one simple pole on the imaginary axis inside 
the physical sheet at 9 = iir/2 + invKi which is associated to a soliton produced as an 
n th -breathers-soliton bound state. 



6. Breather-breather amplitudes 

We now proceed similarly as in the previous section and compute the scattering amplitudes 
of the breathers amongst themselves. As mentioned, this is a very interesting sector as, 
under certain circumstances which will be specified below, it usually closes independently 
from the remaining sectors of the model. Similarly as in (^7l) we exchange now two of the 



auxiliary states (5.1) 



Z n \Gii 92)Z m (9s, 94) — Sb n b rn {9i,92,93,9i)Z m (9s,04)Z n (9i,92) ■ (6.1) 

This is a somewhat cumbersome computation and we will not report here the analogue 
expression of ( |5.3D , since it involves 64 terms, each one of them consisting of a product of 
four amplitudes and four creation operators. As in the previous section, in the next step we 
have to go on-shell by specifying the fusing angles as in (|5,2j). We choose for this purpose 
01 = 00- #2iA 02 = 0o+ 0£l/2, 3 = 0q- 0£/2, 9 4 = 9' + 921/2 with 0% = i(ir-2mvK e ), 
9™ x = i(7r - 2nvK e ), 9 = 9 -9' such that 

013 = + i( n — Tn)vKi, #14 = — iri + i(n + m)vKi, (6-2) 
024 = + i(jn — n)vKi, 023 = 9 + iri — i(n + m)vKg. (6-3) 

For this choice of the fusing angles there are several on-shell identities of the type (5.4)- 
Q5.6j ). Extracting here the common factors of dij we compute 

(-l) m (6 2 4 C13 CM + 614 d 2 3 d-24,) + (-1)™ (6l3 C14 C 2 4 + &14 ^13 ^23) 

+C13 C14 C23 C 24 + 6l4 £>23 + (-l) m+ ™ ( b U 6 2 3 ^13 ^24 + &13 &24 C14 C 23 ) = (6.4) 

(-l) m (6 23 dl3 dl4 + 613 C23 C 24 ) + (-1)™ (^24 C13 C 23 + ^23 ^14 ^24) 

+di3 du d 2 3 d-24 + h3 i>24 + (-l) m+n [d\4 d 2 3 + C13 C24) 
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and 

(-l) m (613 o 24 du + di 3 d 2 3 d 2 4:) + (-1)™ (£13 di4 + <%) 

+&13 C 2 3 C 2 4 du + 2 3 ^13 + (-l) m+n (6 2 3 ^24 + 24 &13 <?23 ^u) = , (6.5) 

(613 614 C 24 + Cu dl 3 J23) + ("!)" ( &14 &24 C13 + C i4 d 2 3 0^4) 

+623 C14 dl3 ^24 + Z>13 &14 &24 C 2 3 + (-l) m+n (S23 C14 + &14 &13 C23 C24) = . (6.6) 

The relations ( |6. 5| ) and ([0]) lead to a cancellation of the backscattering terms in an 
analogous fashion as in ( |5.3[ ). With ( |6,4| ) we indeed end up with a diagonal scattering 
matrix 

Z n (9 )Z m (9' ) = S bnbm (9)Z m (9' )Z n (9 ) (6.7) 

where 



Sb n b m {0) =013014 023 024 C13C14C23C24 + 014^23 + (-1)" m (^14 023^13 ^24+^13 »24 Cl4C 23 



+ (-l) m (6 2 4 C13 Ci4+6i 4 (i23 tfet) + 



1 - fen 2 -sn 2 ( — + (n + m)i^ 



l) n (6 M 0?i 3 d 2 3 + 013 C14 C 24 

1 



(6.8) 



fen 2 — sn 2 



— + (n + m)K e + - 



sn(iO/u — ir/v + (n + m)K() 
sn{i9 /u + ir/v + (n + ra)K^ 



-013014 023 O24 



and 



013014 023 024 = $13$14$23$24 



iWl + f + S + . 



+ 1) 



■p / — m n 
1 ;}2 o o 



oo n— 1 



nn 

k=l 1=1 
oo n— 1 

nn 

k=0 1=1 
oo m— 1 

nn 

k=l 1=0 
oo m— 1 

nn 

k=0 1=0 



III I 7fc 

T 2 



r>(i + 2 

Z-6 



m _|_ n 



2 



+ f)I>( 



-m 
2 



-m 
2 



2 J 
A],}2 



(6.9) 



(6.10) 



2 "i" 2 



Z+0- fcA+AW^ 



2 



+ 1 



1 + + 



2 ^]cj 2 [ ™ 2 + ' 



A; A] ^2 



m i n 
2 "r" 2 



fcALa[ 



m 
" 2 



f + Z + 



^A] £ 



I 2 ~ l ~ 2 



A:A + A]g 2 [-f - 8 + Z + #-fcA]< 



f + §- Z + #- £:A + %[-f-§ + Z 



fc A]g2 



2 2 



Z + 



k A] ^2 



A; A] (.2 



[m _i_ n ; 

I 2 ~ l ~ 2 



| - fcALa[-f - 2 + 1 + 9- |- ifeAU 



The latter expression (6.10) is tailored to make contact to the expressions in the literature 
corresponding to the trigonometric limit. 

The matrix S bnbm (6) a ls° exhibits several types of poles, i) simple and double poles in- 
side the physical sheet beyond the imaginary axis, ii) double poles located on the imaginary 
axis, iii) simple poles in the non-physical sheet and iv) one simple pole on the imaginary 
axis inside the physical sheet at 9 = 9 b = iv{n + m)K^ which is related to the fusing 
process of two breathers b n + b m — ► b n+m . To be really sure that this pole admits such 
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an interpretation, we have to establish that the imaginary part of the residue is strictly 
positive, i.e. 



-i lim (e-e b )s bnbm (e)>o. 

o— >t>b 



(6.11) 



Since the scattering process is parity invariant, we have S bnbm = Sb m 6„, such that we can 
choose without loss of generality n > m. With this choice we compute 



Res w0)=;(-ir +m ^^ 



sinh 



tt— — (n + m) 



fen 4 — ^ ( 5- 2 ( n + m )(™+ m + A + 1 ) 



1 — £sm — l q 
v J 



n-l 



1"^ _j_ YTL Zl ~2 ^ [77, —I— 777, /] ' 

x$ [in] $ [2inuK e ] $ [2invK e ] $ [2i(n + m)zAFQ - ivr] — S S- ] [ - 



z=i 



[n — A]g2 [—ri — A; A] 32 [n + 5 — A; A]g2 [— n + ^ — k X]g2 [n + m — k X]^ 
1 \ [m — k A]^2 [— m — k X] q 2 [m + ^ — k X] q 2 [— m + ^ — k X]g2 [n + m + ^ — k A]^2 

[-n - m - k % [| — A; X}% "^[n + m-l-k X]% [I - n - m - k X}% 



n 



-n 



[— Z - A:A]|j[Z - A; A]? 2 
[^ + ra + m- Z- /c A]|, [f-n-m + Z- A; A] J, \ 



[A_ Z _ A . A] 2 2[ A + ,_ A;A] 2 2 



(6.12) 



The first line in ( 5.12 ) equals i(— l) n+m n with k £ M + . Noting that the functions $ with 
the above arguments are positive real numbers, the second line in Q6.12j ) is {—\) n+m k' with 
n 1 G M + . Recalling finally that 11 + ra < -A/2 we deduce that n&li( ) 6 ^ + sucn tnat 



( |6.11[) is indeed satisfied. 

Due to the factorizability of the theory this fusing process can be associated in the 
usual fashion to a bootstrap equation. For consistency the following equations have to be 
satisfied 



Sib n+m (0) = Sib n (0 + iumKi)Sib m (9-iunK e ) for Z e {b^, s, s} ;k,m + n < n max . (6.13) 



With some algebra we verified (3.13) for the amplitudes derived above Q5.9| ) and ( |6.9| 



7. Reductions of the Z 4 -model 



The elliptic sine-Gordon model can be considered as a master theory, which contains many 
other theories as submodels. By choosing various specific values for the two free parameters 
of the model i and v, one obtains these different types of theories. For the different choices 
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we have the following interrelations 

1 jv — •r2nK^/-K+2imKi _ ( /n 



elliptic sine-Gordon 



elliptic Z^^-affine Toda theory 



m ^ OJ -> 

I 



1 



free theory 



T 

l/u — ► ioo 



m = OJ -> 
I 



(7.1) 



/ 



sine-Gordon 



minimal -D^j^-affine Toda theory 



l/u—*n 

Let us discuss this schematic diagram in more detail: 
7.1 Trigonometric limits 

As was already stated in , when carrying out one of the trigonometric limit I — > for the 
Z4-model for the amplitudes (|4.3| )-( fP| ), one recovers the soliton sector of the sine-Gordon 
model. In our formulation this can be seen directly, as we only have to use the relations 
( |2.9| ) and ( 2.22|) for the Tg- and sn- functions and note that lim^_>o 3> = 1. We have employed 
here the somewhat generalized conventions of the sine-Gordon model formulation in p2| . 
To make contact with the infinite product representation as presented in the literature (see 
equations (2.16) and appendix C in |32[|), we also have to use the identity 



n 

k=0 



r(f + a + 7 )r(f + a 



7) 



n 

fc=0 



r 2 (^ + T) 



r(f + <f + ^)r(f + 



2 i 



(7.2) 



after the limit i — > is performed. As we demonstrated above, when relaxing the constraint 
v > ir/2Ki, the Z4- model exhibits also a breather sector. It is easily seen that the 
corresponding amplitudes can also be obtained in that limit. The expressions for the 
soliton-breather amplitude ([T^) and ( 5.1C| ) are tailored in such a way that we obtain the 
corresponding amplitude in the sine-Gordon model (see equation (20) in |27j) upon the use 
of B , Q and 

sin(x7r/A) -^r (x — kX + A)(— x — kX) 



sin(yvr/A) JJ (y - kX + X)(-y - kX) 



(7.3) 



Similarly, we recover the breather-breather amplitude of the sine-Gordon model (see equa- 
tion (22) in |27f| ) using in addition to (fO| ) also 



n 



(y-fcA -iK-y 



fcA- |) 



cos(x7r/A) 

cos( y7 r/A) ( s _ fe A- D^-fcA-f) 



(7.4) 



There is of course the other trigonometric limit £ —* 1, which one could in principle 
compute by exploiting the relations between the q and g-deformed quantities mentioned in 
section 2. However, just by considering the pole structure ( |4.8| ) and recalling (|2.3| ), we see 
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that there are no poles left inside the physical sheet which could produce a bound state, 
such that this theory will only possess a soliton sector. We will not consider this case here. 

Whereas these limits served essentially only as a consistency check, the next one will 
lead to a new type of theory. 

7.2 Diagonal limit 



For the sine-Gordon model it is well known ]4l|, 42 that in the limit u — ► 1/n the backseat 



tering amplitude vanishes and one obtains a diagonal S-matrix which can be identified with 
a minimal D^\- affine Toda field theory. We observe, that there is an analogue to this 
behaviour in the Z4-model as also in this case the backscattering amplitudes vanish in the 
limit 

l/v -» l/v n ,m = (2nK e + i2mK 1 _ l )/TT . (7.5) 

We find 

lim c(9) = for n,m £ Z, (7-6) 

lim d(9) = for n,m£Z. (7.7) 



For the remaining amplitudes (|4.3| ) and (|4.4| ) in the soliton sector we compute for this limit 

lim b(8) = (-l) n+1 a d {9) for n,m£Z, (7.8) 

lim a {0) = a d {6) forn,meZ. (7.9) 



where 

oo n— 1 



a (e) = TT TT l~U n + mr ) + ( 2fc + l ) n + + mr ) + ( 2fc + l ) n ~ l k 2 (7 w) 

l} \},[i{n + mT) + {2k + l)n + l] e [-t{n + mT) + {2k + l)n-l]^ K ' 

x + |]T 1 h|rr^l - n + (| - m)r + £(n + mr)]r 4 [n + (m - |)r - |(n + mr)] 
S r { [l - n + (I - m)r]r^[n + (m - i)r]r 5 [l + | - ^(n + mr)]r$[-§ + £(n + mr)] 

This is obtained when replacing in 



lim 0— >— — (n + mr) and lim A -> -2(n + mr), (7-11) 



ITT 



and exploiting the property ( |2.15| ) thereafter. Of course one may also carry out the limit 
( |7.5| ) for the soliton-breather and for the breather-breather amplitude. We do not report 
those expressions here as they are quite obvious, unlike in ( [7.10 ), where several cancellations 
could be carried out and the infinite product in the Tg could be turned into products in 
[ ]q2. The scattering matrix obtained in this way belongs to a new type of theory. One 
should mention that the case m ^ 0, which complexifies the coupling constant is most 
likely only of a formal nature, but we expect the case m = to be a meaningful theory. 
However, this requires more analysis especially since there are additional Tachyons present. 

Let us now carry out the limit i — * and verify that the above mentioned diagram 
( |7.1| ) is indeed commutative. For the cases m = 0, n > 3 we reproduce the entire scattering 
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matrix of the minimal D^-aiEne Toda field theory. In particular when carrying out the 



limit in (7.10) we find, upon the use of (ITT 



r , m TT "tt 1 + 2fc + 1) + l][n(x + 2k + 1) - I] 

,J^=o ad{6) = 11 11 [„(x + 2/b + l) + i][n(-x + 2fc + l)-Z] (7 - 12) 

fc=0 t=0 

= n n 1 sinhI(^ + M) 



sinh 



which coincides with the S n +i,n+i = S^-amplitudes in the minimal L'„_j_ 1 -affine Toda 
field theory. For the cases m = 0,n = 1,2 we reproduce the minimal A^,A^- 
affine Toda field theories. We remark that these correspondences hold up to a change of 
statistics, that is some amplitudes are only recovered up to a factor of —1, which changes 
bosonic to fermionic statistics or vice versa. These facts make it natural to call the above 
theories elliptic versions of the associated limiting theory. As we mentioned above, taking 
the breather sector alone constitutes a consistent theory in itself. For q = l,m = we 
have minimal A%^_ 2 C -D^^-affine Toda field theory, a property which transcends also 
into the elliptic version. A final comment is related to the center of the diagram. We find 
lim^i^^o o<i(0) = 1 simply due to the property (|2.7[), 

8. Conclusion 

We demonstrated that when one relaxes the constraint on the coupling constants, one can 
construct a consistent breather sector for the elliptic sine-Gordon model. The scattering of 
the breathers amongst themselves and with the soliton sector satisfies a bootstrap equation 
related to the fusing of two breathers to a third. For the formulation of the scattering 
amplitudes we used as natural objects q-deformed functions. Roughly speaking one replaces 
in the amplitudes in the soliton-antisoliton sector the infinite products of Euler's gamma 
functions by a q-deformed version and in the infinite product in the breather sector integers 
by q-deformed ones. This tailors the models automatically in a form which allows to carry 
out various limits. Instead of carrying out the bootstrap analysis one could alternatively 
take a spin chain as a starting point and use a method based on the algebraic Bethe ansatz, 
pursued for instance in [jSjfl , to compute the breather S-matrix amplitudes. It would be 
interesting to compare that approach with our findings. 

In the diagonal limit we obtain an interesting new theory, which can be viewed as an 
elliptic generalization of the minimal D^^-affine Toda field theory. In [44] we proposed a 



procedure which also lead to elliptic generalizations of theories whose scattering amplitudes 
can be expressed in terms of trigonometric functions. The theories obtained in that fashion 
were, however, of a quite different nature. The procedure in [44] works strictly on the 



principle that there are no redundant poles present in the amplitudes, such that as a 
difference the string of tachyonic states which was encountered here was confined to the 
non-physical sheet where they can be viewed as unstable particles. Further investigation 
is needed to clarify more the interpretation of these models and in particular to establish 
whether they posses a meaningful conformal limit [45|. 
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From a mathematical point of view it will be interesting to generalize the method of 
section 3 to affine Weyl groups of higher rank. 
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